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We  have  performed  an  extensive  theoretical  study  on  die  radiation  properties  of  DQW  structures,  such  as 
gain,  loss,  and  efficiency.  We  have  also  performed  an  extensive  study  using  magnetotunneling 
spectroscopy  to  accurately  determine  the  relative  subband  levels  in  the  DQW  structures. 


14.  SUBJEa  TERMS 

DTI®  QUALITY  INSPECTTED  5 

15.  NUMBER  OF  PAGES 

IS.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 

18.  SECURITY  aASSIFICATION 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 

C* 

20.  LIMITATION  OF  ABSTRACT 

UL 

TOO  7-5191 

UNCLASSIFIED 


FAR-INFRARED  LASERS 
Final  Report 


Tnhie  of  Contents 


Accesion  For 


NTIS  CRA&I  J 

OTIC  TAB 
Unannounced 
Justification _ _ _ 

□ 

By  -  -  — 

Distribution/ 

Availability  Codes 

Dist 

tL 

Avail  a 
Spe 

nd/or 

cial 

Page 


I.  Abstract 

IL  Introduction 

III.  Gain,  loss,  and  lasing  threshold  for  DQW  structures 

III.  1  Radiation  properties  and  gain  of  a  single  DQW  device 

111.2  Scattering  and  lasing  threshold  current  density 

111.3  Magnetic  suppression  of  nonradiative  intersubband  scattering 
in.4  Magneto-tunneling  spectroscopy 

III.5  List  of  participating  scientific  personnel 

s 


2 

3 

4 

11 

14 

16 

20 


I.  Abstract 


This  project  is  to  develop  THz  solid-state  lasers  using  multiple-quantum-well  (MQW) 
structures.  Based  on  work  in  our  group  and  other  groups  outside  MIT  in  the  last  few  years,  we  have 
focused  our  effort  on  using  double  quantum-well  (DQW)  structures  to  achieve  the  goal.  In  essence,  a 
DQW  device  is  a  voltage-tunable  two-level  system,  as  far  as  the  interwell  intersubband  transition  is 
concerned.  Thus,  compared  to  intrawell  intersubband-transition  devices  that  have  been  proposed  and 
studied,  the  DQW  device  has  a  greater  tolerance  for  design  and  fabrication  errors.  This  is  a  major 
advantage  over  Bloch  oscillators  using  superlattices,  in  which  the  unavoidable  formation  of  high-field 
domains  will  lead  to  a  significant  misalignment  of  energy  levels  through  the  superlattice.  Since  the 
relative  energy  difference  between  the  lowest  two  levels  in  a  DQW  structure  is  voltage  tunable,  the 
same  structure  can  also  be  used  for  far-infrared  detection  as  well  as  generation,  if  the  initial  state  lies 
below  the  final  state.  In  the  current  three-year  project,  we  have  explored  the  potential  of  far-infrared 
generation  and  detection  using  DQWs.  We  have  performed  an  extensive  theoretical  study  on  the 
radiation  properties  of  DQW  structures,  such  as  gain,  loss,  and  efficiency.  We  have  also  performed  an 
extensive  study  using  magnetotunneling  spectroscopy  to  accurately  determine  the  relative  subband 
levels  in  the  DQW  structures. 


s 
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II.  INTRODUCTION 


Far-infrared  or  THz  frequencies  are  among  the  most  underdeveloped  electromagnetic  spectra, 
even  though  their  potential  application  in  spectroscopy,  astrophysics,  plasma  diagnostics,  remote 
sensing  and  imaging,  and  space-based  radar  and  communications  are  obviously  great.  This 
underdevelopment  is  because  the  far-infrared  or  THz  frequency  range  falls  between  two  other 
frequency  ranges  in  which  conventional  semiconductor  devices  have  been  successfully  employed.  One 
is  the  microwave  and  millimeter-wave  frequency  range,  and  the  other  is  the  near-infrared  and  optical 
frequency  range.  Semiconductor  devices  that  utilize  the  classical  diffusive  transport  of  electrons,  such 
as  diodes  and  transistors,  have  a  high  frequency  limitation  of  several  hundred  GHz.  This  is  limited  by 
electron  mobility  and  the  smallest  feature  size  that  can  be  fabricated.  Semiconductor  devices  based  on 
quantum  mechanical  interband  transitions,  however,  are  limited  to  frequencies  higher  than  those 
corresponding  to  the  semiconductor  energy  gap,  which  is  higher  than  10  THz  for  most  bulk 
semiconductors.  Therefore,  a  large  gap  exists  from  100  GHz  to  10  THz,  in  which  very  few  solid-state 
devices  are  available. 

Semiconductor  quantum  wells  are  human-made  quantum  mechanieal  systems  in  which  the 
energy  levels  can  be  chosen  by  changing  the  dimensions  of  the  quantum  wells.  Long-wavelength 
photoelectric  devices  that  utilize  the  intersubband  transitions,  such  as  far-infrared  lasers,  have  been 
proposed.  Significant  progress  was  made  recently  toward  this  goal.  Large  oscillator  strengths  of 
intersubband  transitions  have  been  observed  in  far-infrared  absorption  speetroscopy  experiments.’^’* 
An  intersubband  spontaneous  emission  with  a  power  level  of  ~10‘^  W  has  been  observed.^’^®  Despite 
the  encouraging  preliminary  results,  far-infrared  lasing  has  not  yet  been  achieved,  even  though 
unipolar  lasers  based  on  intersubband  transitions  have  been  developed  at  near-infrared  frequencies  (A, 
«  4  ^im).ll  Two  main  obstacles  are  especially  difficult  to  overcome  at  far-infrared  frequeneies.  The 
first  one  is  the  extremely  unfavorable  ratio  of  the  radiative  relaxation  rate  1/Xspon  over  that  of  the 
nonradiative  rate  1/Xac  +  1/xlo  +  which  is  due  to  acoustic-phonon,  LO-phonon,  and  electron- 

electron  scattering.  Specifically,  for  intersubband  transitions  within  a  single  quantum  well  (intrawell)  at 
THz  frequencies,  Xspon  ~  1  fis,  Xac  -  100  ps,  Xlo  <  1  ps.  and  Xg-e  ~  1  ps  for  a  earner  density  of 
10l2/cm2  For  interwell  transitions,  these  rates  will  be  Educed  by  the  tunneling  probability  T  through 
the  middle  barrier;  except  e-e  scattering  rate  is  reduced  by  T^,  since  two  electrons  must  traverse  the 
barrier.  The  seeond  obstacle  is  the  mode  confinement.  The  dielectric  waveguide  technology  that  was 
successfully  used  at  near-infrared  frequencies  has  a  low  loss  of  less  than  10  cm"^.  While  the  plasma 
confinement^  used  at  far-infrared  frequencies  has  a  heavy  loss  of  several  hundred  cm  due  to  the  low 
mobility  associated  with  heavy  doping  in  the  plasma  region. 
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In  the  current  project,  we  have  explored  to  use  many  identical  DQW  structures  that  are 
connected  in  series  in  order  to  overcome  these  two  difficulties.  By  using  such  structures,  the  required 
inverted  population  density  can  be  spread  out  to  many  DQW  modules  in  order  to  keep  the  two- 
dimensional  density  low  (for  example,  l.SxlO^Vcm^)  in  each  module.  In  this  way,  the  difference 
between  the  Fermi  energy  and  the  lower  subband  level  will  be  below  ficoLO*  so  that  the  LO-phonon 
scattering  will  not  occur.  Furthermore,  the  electron-electron  scattering  will  be  reduced  at  low  carrier 
densities.  We  have  also  applied  a  longitudinal  magnetic  field  (along  the  current)  to  further  suppress  the 
nonradiative  relaxation  processes  in  order  to  improve  the  radiation  efficiency.  ^2. 13  xhe  active  region's 
greater  thickness  of  many  DQW  structures  will  enhance  the  mode-confinement  factor,  thus  reducing 
relative  loss  due  to  the  plasma  walls. 

III.  Gain,  loss,  and  lasing  threshold  for  DQW  structures 

III.l  Radiation  properties  and  gain  of  a  single  DQW  device^'* 

The  schematic  of  a  single  DQW  structure  is  shown  in  Fig.  1(a).  The  conduction  band  profile, 
subband  structures,  and  carrier  densities  are  calculated  self-consistently  from  Schrbdinger  and  Poisson 
equations,  after  typical  five  iterations.  Fig.  1(b)  shows  the  transmission  coefficient  through  such  a 
structure,  whose  peaks  correspond  to  the  subband  levels  in  Fig.  1(a).  Clearly,  the  intersubband  energy 
difference  E2  -  Ei  (corresponding  to  the  radiation  frequency)  can  be  tuned  by  a  bias  voltage  between 
the  emitter  and  collector.  The  most  important  parameter  that  determines  the  radiation  properties  (due  to 
the  l—>2  intersubband  transition)  of  this  DQW  structure  is  the  dipole  moment  (llzl2),  or  equivalently, 
the  oscillator  strength  fi2,  which  is  the  ratio  of  the  quantum  mechanical  radiative  decay  rate  over  that  of 
a  classical  Hertzian  dipole  at  the  same  frequency. 

Before  proceeding  with  numerical  calculations,  it  is  instructive  to  gain  insight  into  the 
functional  relationship  of  the  oscillator  strength  and  the  energy  splitting  for  a  given  anticrossing  gap  Aq 
(which  is  the  minimum  value  of  IE2-E1I).  This  is  done  by  analytically  solving  the  problem  using  a 
coupled  two-level  model.  Express  the  stationary  eigenstates  II)  and  12)  in  terms  of  superpositions  of 
states  localized  in  the  left  and  right  well,  IL)  and  IR), 

ll)  =  CLML)-hCRMR),  s 

and  12)  =  cl^  IL)  +  cr^  IR) . 

Here,  Icl.r^’^I^  is  the  probability  of  finding  a  particle  in  well  L  (or  R)  for  the  eigenstate  II)  (or  12)).  The 
resulting  Heisenberg  equation  becomes, 

1,2 

r  E.  -  El, 2  -  eV/2  -Ay2  T  f 

L  Ao/2  Eo-Ei,2 +eV/2  J  L 
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Figure  1.  (a)  Conduction  band  profile  and  the  square  of  wavefunctions  at  a  bias  voltage  of  240  mV  for  a  typical 

DQW  structure.  The  wider  emitter  well  (on  the  left)  is  to  ensure  that  a  large  carrier  density  can  be  achieved  at  a  desired 
intersubband  spacing  Ej  -  Ei.  (b)  Transmission  coefficient  T(E)  as  a  function  of  the  electron  incident  energy.  The  peak 
positions  in  T(E)  correspond  to  the  subband  levels  in  (a).  The  resonance  widths  in  T(E)  dctennine  the  escape  lifeumes. 

Here,  E©  is  the  eigenenergy  in  well  IL>  when  it  is  isolated  from  the  well  IR),  V  is  the  potential 
difference  between  the  two  wells  at  the  well  centers,  '^e  solution  of  Eq.  (1)  is 

E2.1  =  Eo  ±  V  (Ao/2)2  +  (eV/2)2  ,  (2) 

and  the  field-induced  energy  splitting  A  =  IE2  -  Eil  =  V(Ao)2  +  (eV)2  .  A  as  a  function  of  the  bias 
voltage  V  is  shown  in  Fig.  2.  The  energy  splitting  takes  the  minimum  value  Aq  at  the  anticrossing 
point,  above  which  the  energy  of  the  state  that  is  mainly  localized  in  the  left  well  becomes  higher  than 
that  of 
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Figure  2.  (a)  Energy  level  diagram  as  a  function  of  applied  voltage,  which  exhibits  an  anti-crossing  with  a  gap 

Aq.  At  low  voltages  the  device  can  be  used  as  a  tunable  detector,  as  indicated  in  (b),  and  beyond  die  anti-aossing  it 

functions  as  a  field-tunable  emitter. 


the  state  that  is  mainly  localized  in  the  right  well,  shown  in  Fig.  1(a).  Since  the  state  mainly  localized  in 
the  left  well  couples  strongly  to  the  emitter  (which  always  has  a  higher  Fermi  energy  than  the 
collector),  and  is  therefore  the  initial  state  that  is  filled  with  electrons.  Fig.  2  indicates  that  the  relative 
energy  difference  between  the  initial  and  final  states  (not  only  the  amplitude  but  also  the  sign)  can  be 
tuned  by  the  bias  voltage.  Consequently,  the  same  DQW  structure  can  be  used  as  both  a  detector  and 
an  emitter.  As  a  detector,  the  initial  state  lies  below  the  final  one;  while  as  an  emitter,  the  initial  state 
lies  above  the  final  one. 

Semiquantitatively,  the  coupling  strength  of  the  DQW  to  the  emitter  and  collector  reservoirs  can 
be  estimated  from  the  probability  of  finding  a  particle  in  the  left  (emitter)  and  right  (collector)  well,  and 


they  are. 


and 

where  9  is  defined  by 


lcj^l,2|2 

lcj^l,2|2 


1  -b  e±0  ’ 

1  -  IcL^’^P  , 


sinh(0/2)  =  eV/Ao . 

These  probabilities  are  plotted  in  Fig.  3.  At  the  anticrossing  gap  A 


Aq.  Icl*’"I^  =  IcR^^P  =  0.5,  the 


wave  functions  for  states  II)  and  12)  are  completely  delocalized.  The  two  wells  have  the  strongest 


coupling  at  this  point.  As  A  deviates  from  the  anticrossing  gap  Aq,  however,  the  probability  deviates 
from  0.5  and  the  wave  functions  are  localized  in  either  the  left  or  the  right  well.  This  is  a  lield-induced 


Figure  3.  Probabilities  Icl^  as  functions  of  energy  splitting  A  and  anti-crossing  gap  Aq. 

Wannier-Stark  localization.  As  shown  in  Fig.  3,  the  probability  lc|2  becomes  approximately  90%  (or 
10%)  at  a  normalized  A/Ao  =  2.  Thus,  the  decoupling  of  the  two  wells  is  a  rather  rapid  function  of  the 
energy-splitting  A/Aq.  This  decoupling  has  the  benefit  of  making  selective  population  and  depopulation 
easier  (since  the  upper-level  state  is  mainly  coupled  to  the  emitter  while  the  lower-level  state  is  to  the 
collector).  However,  it  also  has  the  adverse  effect  of  reducing  the  dipole  moment  (or  the  oscillator 
strength)  of  the  I1)-^I2)  intersubband  transition.  Therefore,  we  usually  consider  the  case  of  A/Aq  =  2-4 
in  our  design  analysis  for  the  best  trade-off  between  selective  injection  and  a  large  oscillator  strength. 
Analytically,  the  dipole  moment  can  be  shown  to  be 


I  (1  1  z  I  2  )  |2 


W2Ao2 
4  ^2  ’ 


(3) 


where  W  is  the  distance  between  the  centers  of  the  f^o  wells.  Consequently,  the  oscillator  strength 


becomes, 

fl2  =  ^IE2-EiI|(11zI2)12 

_  2m*  W2  Ao2  (4) 

~  fi2  4  ^  • 
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Fig.  4  shows  the  oscillator  strengths  as  functions  of  the  emission  frequency  A  =  IE2  -  Eil  (1  THz  =  4 
meV),  with  the  anti-crossing  gap  Aq  as  a  parameter.  Fig.  4  also  shows  the  required  barrier  width  and 
the  A1  concentration  for  a  chosen  value  of  Aq  for  GaAs/AlGaAs  quantum  wells  through  the 


relation,^5,l6 


where 

d  and  Eb  are  the  barrier  width 


27t 


—  sk  ^ 


_  ^-Indfk 


X  =  h/V2m*EB  , 
and  height,  respectively. 


Figure  4.  (a)  Contour  plot  of  tlie  antiaossing  gap  Aq  as  a  function  of  tlie  barrier  height  and  width.  The  right  axis 

presents  the  A1  concentration  as  a  function  of  the  barrier  height  in  tlie  GaAs/AlGaAs  material  systems,  (b)  The 
oscillation  strength  as  a  function  of  emitted  photon  energy  A  and  witli  die  anticrossing  gap  Ay  as  a  parameter. 
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To  ensure  that  the  greatest  energy  difference  between  electrons  in  subband  12)  and  II)  is  below 
the  LO  phonon  energy  of  fiWLO  =  36  meV,  we  have  chosen  our  frequency  to  be  about  4  THz,  which 
corresponds  to  16  meV  intersubband  spacing.  In  this  way,  we  can  safely  populate  the  second  subband 
up  to  Ef  -  E2  =  15  meV  (corresponding  to  a  5x10^  Vcm^  carrier  density)  without  opening  up  the  fast 
LO-phonon  emission  channel.  Given  A  =  15  meV,  and  if  we  choose  Aq  to  be  4.2  meV,  Fig.  4  shows 
an  oscillator  strength  close  to  0.20,  which  is  an  acceptable  value.  The  contour  plots  in  Fig.  4  indicate 
that  a  barrier  with  approximately  25%  A1  concentration  and  a  width  of  24  A  will  provide  the  desired 
Aq.  These  material  parameters  can  be  controlled  in  a  MBE  growth  process. 


5  10  15  5  10  15 

Energy  A  |meV]  Energy  a  [me V) 


Figure  5.  Spontaneous  emission  lifetimes  for  a  3D  and  2D  optical  mode  density  for  interwell  (expressed  as 

symbols)  and  intxawell  (solid  lines)  transitions  as  a  function  of  emission  frequency  and  with  as  a  parameter. 
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From  the  oscillator  strength,  it  is  straightforward  to  calculate  the  intersubband  spontaneous 


radiative  lifetime  tspon 


yi  T  2  ^  1 

4  e  L  c  1  determines  the  gain  of  the  lasing  media.  Here 

fl2 


we 


0) 


assume  the  radiation  is  confined  between  two  parallel  plasma  walls  with  a  separation  of  L.  Assume  L 
=  3  ^tm,  Tspon  is  plotted  in  Fig.  5  as  a  function  of  the  energy  splitting  with  Aq  as  a  parameter.  In 
comparison,  the  radiative  lifetime  of  a  Hertzian  dipole  emitting  in  2D  is  also  plotted  in  Fig.  5,  which 
shows  a  l/co  dependence  (in  contrast  with  l/co^  dependence  of  a  3D  dipole).  The  radiative  lifetime 
Tspon  is  almost  independent  of  frequency  because  of  the  l/m  dependence  of  the  oscillator  strength  (see 
Eq.  (4)).  For  Aq  =  4.2  meV,  this  lifetime  is  approximately  6  [is. 

The  gain  of  a  medium  with  an  inverted  population  An  can  be  calculated  from  the  well-known 
Einstein  relation  between  spontaneous  and  stimulated  transition  rates,  and  is  given  by, 

^  _ Jz _  (5) 

^  hot  tac  dtrerTspon  1  +  (O)  - 

where  tac  is  the  thickness  of  the  active  region  with  an  inverted  population,  ttot  is  the  distance  of  photon 
confinement,  cDq  =  (E2  -  Ei)/fi  is  the  angular  frequency  at  resonance,  tzT2  is  the  inverse  of  the  full 
linewidth  of  the  spontaneous  emission,!^  and  Er  is  the  relative  dielectric  constant  of  the  quantum  well 
material.  At  the  lasing  threshold,  the  gain  equals  the  total  loss,  which  is  mainly  due  to  free-carrier 
absorptions  in  the  lossy  plasma  walls.^  This  loss  has  been  calculated  before  and  is  given  by. 


a  = 


1  1 


hot  Sj 


(OpT 


hot’ 


(6) 


where  Sz  and  Sx  are  the  z-  (along  the  growth  direction)  and  x-component  of  the  Poynting  vector, 
respectively;  cOp  is  the  plasma  frequency,  and  t  is  the  scattering  time  in  the  plasma.  Equating  Eqs.  (5) 
and  (6),  we  obtained  the  lasing  threshold  density  of  population  inversion,  which  is  plotted  in  Fig.  6  as 
a  function  of  emission  frequency  for  two  different  doping  levels  in  the  plasma  region.  Again,  the  anti¬ 
crossing  gap  is  a  parameter  in  these  calculations.  Clearly,  there  is  some  benefit  in  using  heavily  doped 
plasmas,  whose  shorter  penetration  depths  (corresponding  to  tighter  mode  confinements)  more  than 
compensate  for  the  shorter  scattering  time  t  at  a  higher  doping  level.  The  approximate  0)^2  increase  of 
the  threshold  density  is  due  to  the  increase  of  th^))lasma  loss,  which  is  analogous  to  the  surface 
resistance  of  metals. 

If  we  use  DQW  structures  with  a  Aq  =  4.2  meV  and  a  doping  concentration  of  2xl0l8/cm3  in 
the  plasma  regions,  the  threshold  density  will  be  approximately  lO^^/cm^.  If  all  these  electrons  are 
placed  in  a  single  DQW  device,  the  Feimi  energy  will  be  approximately  30  meV,  which  will  allow  fast 
LO-phonon  scattering  to  occur.  Furthermore,  the  electron-electron  scattering  time  will  be  shorter  than  1 
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ps.  These  two  extremely  fast  relaxation  processes  will  make  the  required  inverted  population  density 
impossible  to  achieve.  However,  if  we  use  many  identical  DQW  modules  to  spread  out  the  required 
inverted  population.  An  will  be  much  less  than  I0i3/cm2  in  each  of  them.  Consequently,  the  LO- 
phonon  emission  process  can  be  avoided  and  the  electron-electron  scattering  will  be  substantially 
reduced. 


Figure  6.  Threshold  inverted  population  density  required  to  overcome  plasma  confinement  losses  for  two  plasma 

concentrations:  1x10**  cm'^  (a)  and  2x10**  cm'^  (b).  For  these  calculations,  the  linewidth  1/^X2  is  taken  to  be  20 

cm'*  from  Ref.  [10]. 

III. 2  Scattering  and  lasing  threshold  current  density 
The  lasing  threshold  cunent  density  is  simply  given  by 

Jth  =  e(Aniii/N) /iioi , 
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where  N  is  the  number  of  DQW  structures,  and  1/Xiot  =  +  l/'^e-e  is  the  total  relaxation 

rate.  (We  ignore  the  extremely  slow  tspon  radiative  relaxation  here.)  In  the  current  project,  we  have 
made  major  efforts  to  calculate  the  rates  of  all  nonradiative  relaxation  processes.^^  We  summarize  the 
results  in  the  following.  Fig.  7  shows  xlo  as  a  function  of  the  radiation  frequency,  assuming  the 
electrons  in  the  second  subband  have  sufficient  kinetic  energy  to  emit  an  LO  phonon  in  the  case  of  E2  ■ 
El  <  ficoLO-  For  ^0  =  4.2  meV,  xlo  ~  1  PS  for  the  frequency  range  of  interest.  This  extremely  fast 
relaxation  makes  the  radiation  efficiency  low  (-10-7)  and  the  threshold  current  density  prohibitively 
high.  Thus,  this  channel  should  be  avoided  in  any  design  considerations. 


Energy  Splitting  A  [meV] 


Figure  7.  The  GaAs-like  LO-phonon  scaitering  time  Tlo  for  a  DQW  structure  as  a  function  of  subband  spacing  A 

and  at  different  anticrossing  gaps  A^.  The  plot  on  die  right  is  the  radiative  efficiency  XLofXspon  a  function  of  A. 
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Acoustic-phonon  scattering  is  a  much  slower  process  because  of  the  absence  of  singularities  in 
the  density  of  phonon  states.  Fig.  8  plots  Xac  as  a  function  of  frequency.  Again,  for  Aq  =  4.2  meV,  Xac 
~  1  ns.  This  process  is  sufficiently  slow  so  that  it  will  not  require  a  demanding  current  density  to  reach 

the  threshold. 


Figure  8.  The  acoustic-phonon  scattering  time  "tac  as  a  function  of  A  with  Aq  as  a  parameter.  The  plot  on  the 

right  is  the  radiative  efficiency  Tac/Xspon  4ue  to  acoustic-phonon  scattering. 
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Electron-electron  scattering,  which  has  been  largely  ignored  in  most  of  the  previous  work,  was 
found  to  be  the  dominant  nonradiative  relaxation  process  in  this  project.  Although  e-e  scattering  is  an 
elastic  process,  as  far  as  the  total  energy  of  the  two  electrons  is  concerned,  it  can  result  in  a  large 
change  of  moment  for  each  electron.  Therefore,  two  electrons  can  scatter  each  other  from  the  second 
subband  to  the  first  one,  gaining  a  greater  transverse  (in-plane)  momentum  at  the  expense  of 
longitudinal  energies.  These  electrons  will  later  emit  phonons  to  lower  their  energies  to  the  bottom  of 
the  first  subband.  As  far  as  the  radiation  property  is  concerned,  this  scattering  is  equivalent  to  other 
nonradiative  processes,  such  as  acoustic-phonon  and  LO-phonon  emissions,  which  are  inelastic 
processes.  The  estimate  of  the  e-e  scattering  rate  is  by  no  means  trivial  and  some  calculations  remain 
controversial.  Using  the  simplest  Hartree  formalism  and  ignoring  the  exchange  interaction,  we 
estimate  the  e-e  scattering  time  to  be  approximately  30  ps  for  Aq  =  4.2  meV  and  A  =  10  meV,  for  a 
modest  carrier  concentration  of  1.4x101  Vcm2.14  Clearly,  e-e  scattering  is  the  dominant  nonradiative 

relaxation  process  in  the  frequency  range  of  interest. 

For  the  required  threshold  density  of  Anth  =  lOl^/cm^,  if  we  use  70  DQW  modules  to  reduce 
An  to  1.4xl0ll/cm2  for  each  one,  the  lasing  threshold  current  density  from  Eq.  (7)  will  be 
approximately  750  A/cm^.  This  is  feasible,  provided  the  property  of  these  DQWs  is  sufficiently 
uniform  to  evenly  spread  out  the  inverted  population. 

III. 3  Magnetic  suppression  of  nonradiative  intersubband  scattering 

Clearly,  if  we  were  able  to  reduce  the  nonradiative  intersubband  scattering  rate,  especially  the 
LO-phonon  and  electron-electron  scattering,  we  could  reduce  the  lasing  threshold  current  density  Jth 
significantly;  or  to  keep  Jui  constant,  we  could  increase  the  inverted  population  density  significantly, 
and  therefore  using  fewer  DQW  modules  and  still  have  sufficient  lasing  gain.  The  nonradiative 
relaxation  processes  are  so  fast  because  the  subbands  in  a  quantum  well  have  continuum  states  for  all 
the  possible  transverse  momentum  values  (perpendicular  to  the  growth  direction),  instead  of  discrete 
energy  levels.  Thus,  the  conservation  of  energy  and  transverse  momentum  can  be  simultaneously 
satisfied  for  these  processes.  In  a  quantum  dot  with^  a  zero-dimensional  density  of  states,  these 
processes  will  be  significantly  suppressed,  since  the-LO-phonon  process  will  be  prohibited  if  the 
intersubband  spacing  is  different  from  ficoLO-  Electron-electron  scattering  will  also  be  suppressed,  if 
the  lower  and  upper  subbands  have  different  spacings,  as  illustrated  in  Fig.  9.  In  addition,  the 
acoustic-phonon  process  will  be  reduced,  because  the  number  of  available  phonon  modes  is 
significantly  reduced.  However,  a  quantum  dot  contains  few  electrons,  thus  it  will  not  generate  much 
radiation  power.  A  quantum-dot  array  will  inevitably  introduce  inhomogeneous  broadening.  In 
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addition,  an  electiic  contact  to  such  arrays  will  be  diflicult,  thus  making  electric  pumping  a  less  feasible 
approach. 


Figure  9.  (a)  Elecuon-electron  intersubband  scattering  in  a  quantum  well  that  simultaneously  satisfies  the  energy 

and  momentum  conservation,  (b)  Suppression  of  the  e-e  scattering  in  a  zero-dimensional  system. 

Similar  to  the  concept  employed  to  confine  plasmas  in  a  tokamak,  several  groups  have 
proposed  the  use  of  magnetic  fields  to  reduce  the  dimensionality  of  the  electron  system  in  order  to 
suppress  nonradiative  relaxation  processes.^^’^^  Under  a  strong  field  along  the  growth  direction,  the 
transverse  motion  of  electrons  will  be  confined  in  cyclotron  orbits  and  their  energy  spectrum  will  be 
quantized  at  Landau  levels.  This  magnetic  field-induced  transverse  confinement,  along  with  the 
longitudinal  confinement  provided  by  the  quantum-well  structures,  will  effectively  reduce  the 
dimensionality  of  the  electron  system  to  that  of  a  zero-dimensional  system,  similar  to  that  of  quantum- 
dot  systems.  Consequently,  the  detrimental  nonradiative  relaxation  processes  will  be  significantly 
reduced.  In  this  scheme,  the  applied  longitudinal  magnetic  field  will  not  affect  transport  in  that 
direction  or  the  dipole  moment  associated  with  the  transition  between  two  subbands.  Therefore,  the 
field  will  not  have  any  adverse  effect  on  the  radiative  process,  and  will  substantially  increase  the 

radiation  efficiency  ii  =  'Cnonrd/'^spon-  \ 

In  the  current  project,  we  have  designed  and  fabricated  appropriate  MQW  structures  in  which 
the  current  is  limited  by  relaxation  between  two  subbands,  which  is  mostly  nonradiative  below  the 
lasing  threshold.  Under  a  longitudinal  magnetic  field,  this  relaxation  process  and  the  current  will  be 
reduced.  This  will  provide  definitive  evidence  for  the  predicted  magnetic  suppression.  Our  preliminary 
measurements  have  yielded,  at  best,  ambiguous  results  to  support  this  mechanism.!^  Other  groups 
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pursuing  this  goal  have  experienced  similar  results.  We  speculate  that  the  failure  to  observe  the 
magnetic  suppression  is  due  to  broadened  Landau  levels  that  resulted  from  impurity  and  interface 
scattering.  These  broadened  levels  provide  the  states  for  scattering  which  simultaneously  satisfy  the 
energy  and  momentum  conservations.  We  intend  to  improve  this  measurement  by  using  better  MBE 
samples  and  stronger  magnetic  fields  (previous  work  was  conducted  up  to  15  Tesla),  at  which  the 
Landau  gap  is  opened  up  faster  (scaled  linearly  with  H)  than  the  level  broadening  (scaled  as  Hl/2 
[13]).  This  will  reduce  the  number  of  states  for  scattering  to  occur.  At  MIT,  we  have  the  benefit  of  the 
Francis  Bitter  National  Magnet  Laboratory  (FBNML),  which  has  the  high  magnetic  fields  required  for 
this  investigation. 

It  will  also  be  interesting  to  study  the  far-infrared  emission  properties  of  these  DQW  devices 
under  a  strong  longitudinal  magnetic  field.  If  magnetic  suppression  occurs  and  does  not  affect  the 
dipole  intersubband  transitions,  the  emission  linewidth  should  be  narrowed,  due  to  reduced  scattering. 
For  the  same  reason,  the  intensity  at  the  peak  frequency  of  the  emission  spectrum  should  be  increased, 
because  more  photons  will  be  emitted  at  the  same  frequency.  If  established  experimentally,  this  effect 
will  be  a  beautiful  example  of  magnetic  field-induced  reduced  dimensionality,  which  will  play  an 
important  role  in  device  applications. 

III. 4  Magnetotunneling  spectroscopy 

Motivated  by  the  magnetic  suppression  effect,  we  have  performed  extensive  tunneling  studies 
of  DQWs  and  TQWs  (triple  quantum  wells)  in  strong  magnetic  fields  with  both  a  longitudinal  and 
transverse  component.  This  study,  called  inagnetotunnelhrg  spectroscopy,  is  a  very  powerful 
spectroscopic  tool  to  measure  the  subband  levels  and  in-plane  energy-momentum  dispersion  (crucial 
information  for  any  intersubband  devices)  with  much  greater  accuracy  than  a  dc  tunneling  study  alone 
would  yield. 

At  far-infrared  frequencies  of  several  THz,  the  corresponding  intersubband  spacing  is  only  on 
the  order  of  10  meV  (1  THz  =  4  meV).  For  such  a  narrow  intersubband  spacing,  the  dc  tunneling  I-V 
characteristic  measurement  usually  exhibits  many  rich  features  that  will  be  difficult  to  unambiguously 
identify.  Therefore,  it  will  be  difficult  to  determine  the  relative  subband  positions  from  the  dc  tunneling 
I-V  curves  alone.  The  magnetotunneling  spectroscopy  tool  developed  in  this  work  can  solve  this 
problem  easily.^®  The  current  at  a  fixed  bias  voltage  of  MQW  structures  exhibits  an  oscillatory 
dependence  on  the  longitudinal  magnetic  field.  This  magneto-oscillation  is  periodic  with  the  inverse  of 
the  magnetic  field,  as  shown  in  Fig.  10,  and  is  reminiscent  of  de  Haas-Shubinikov  (dHS)  oscillations. 
Similar  to  using  the  dHS  effect  to  measure  Fermi  energies,  we  can  measure  the  intersubband  spacing 
from  the  oscillation  period  with  sub-meV  accuracy  [20]. 
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The  physical  origin  of  this  magneto-oscillation  is  scattering-mediated  inter-Landau-level 
tunneling.  Suppose  the  energy  in  the  first  well  is  Ei  -h  (m  +  l/2)fieB/m*.  and  the  energy  in  the  second 
well  is  E2  +  (n2  +  l/2)fieB/m*,  where  Ei  and  E2  are  the  energies  of  the  subbands  in  the  first  and 
second  well.  In  the  absence  of  LO-phonon  emissions,  energy  conservation  requires  that  the  two 
energies  be  equal.  Thus,  (n2  -  ni)fieB/m*  =  Ei  -  E2,  or  1/B  =  (n2  -  ni)fie/m*(Ei  -  E2).  This  implies 
that  whenever  1/B  is  an  integer  times  the  inverse  of  a  characteristic  field  I/Bq  =  fie/m*(Ei  -  E2),  the 
energy  conservation  is  satisfied  and  the  tunneling  current  will  be  enhanced.  This  gives  the  oscillatory 
behavior  of  the  tunneling  current  as  a  function  of  1/B,  as  shown  in  Fig.  10,  whose  Fourier  transform 
will  have  a  peak  centered  at  the  characteristic  field  Bq. 


Magnetic  Field  B„m  Fundamental  F.eld  0,  m 


Figure  10.  Magneiotunneling  oscillations  (left)  and  their  Fourier  transform  (right)  in  the  current  of  a  DQW  device 

with  a  bias  voltage  from  489  mV  to  524  mV.  The  inset  on  Uie  left  shows  tlie  frequency  of  tlie  magneto-oscillations  as  a 
function  of  the  angle  between  the  magnetic  field  and  current  flow  The  cost))  dependence  is  an  indication  that  the 
magneto-oscillation  is  due  to  2D-2D  tunneling  between  two  wells. 
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Since  the  energy  difference  between  the  two  subbands  is  related  to  Bq  through  Ei  -  E2  = 
lieBo/m*,  magnetotunneling  spectroscopy  provides  us  with  a  very  accurate  method  to  measure  the 
relative  positions  of  the  subbands.  Fig.  1 1  shows  the  subband  energy  levels  as 
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Figure  11.  I-V  characteristic  (a)  and  Uie  subband  energie.^(b)  as  functions  of  the  bias  voltage  at  4.2  K.  The  quasi- 
Fenmi  level  in  tlie  emitter  well  is  used  as  tlie  energy  reference  (i.€’‘£FE  =  0).  The  opcn-circlc  data  in  the  negative  dynamic 
region  (NDR)  was  obtained  by  applying  an  in-plane  field  to  broaden  die  resonance  peak. 


functions  of  bias  voltage  for  a  DQW  device.  These  measured  energy  levels  are  within  5%  from  the 
designed  values,  which  is  remarkable  considering  the  relatively  small  energy  scales  (~10  meV)  for  far- 
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infrared  applications.  This  accurate  determination  of  subband  energy  levels  is  crucial  to  successfully 
develop  far-infrared  sources  with  specific  frequencies. 

Under  a  transverse  or  in-plane  magnetic  field  (with  the  field  perpendicular  to  the  current-flow 
direction),  as  well  as  a  longitudinal  one,  electrons  will  acquire  momentum  in  the  transverse  direction 
while  traversing  the  MQW  structures.  Thus,  the  conservation  of  transverse  momentum  is  no  longer 
valid.  Consequently,  the  applied  transverse  field  broadens  the  resonant  peaks  in  the  tunneling  I-V 
characteristics,  and  shifts  the  position  of  the  resonances  to  higher  bias  voltage,  as  shown  in  Fig.  12. 
Using  this  method,  we  can  perform  spectroscopy  in  the  I-V  curve  region,  where  it  has  a  negative 
dynamic  resistance  (NDR)  that  otherwise  would  be  inaccessible,  due  to  the  instability  associated  with 
the  NDR.  This  information  is  crucial  for  laser  operation,  since  the  NDR  region  corresponds  to  the 
initial  state  with  a  higher  energy  than  the  final  state,  which  is  the  normal  operating  condition  for  lasers. 


kpoVAj 

Figure  12.  Left  plot:  I-V  tunneling  characteristics  near  resonance  of  an  InGaAs/InAlAs  DQW  device  with  an  in¬ 
plane  field  varying  from  0  T  to  12.9  T  in  steps  of  1.17  T.  Right  plot:  In-plane  energy-momentum  dispersion  relation 
whose  curvature  gives  tlie  effective  mass  m*,  which  is  very  close  to  the  effective  mass  of  bulk  InGaAs  materials. 
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In  addition,  by  studying  the  shift  of  the  current  peak  position  as  a  function  of  the  transverse  field,  we 
can  accurately  measure  the  in-plane  energy-momentum  dispersion  relation,  as  shown  in  Fig.  12,  which 
is  another  piece  of  useful  information  for  the  MQW  structures  [18]. 
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